Crucial role of sidewalls in velocity distributions in quasi-2D granular gases 
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Our experiments and three-dimensional molecular dynamics simulations of particles confined to 
a vertical monolayer by closely spaced frictional walls (sidewalls) yield velocity distributions with 
non-Gaussian tails and a peak near zero velocity. Simulations with frictionless sidewalls are not 
peaked. Thus interactions between particles and their container are an important determinant of 
the shape of the distribution and should be considered when evaluating experiments on a constrained 
monolayer of particles. 
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Granular materials can mimic the behavior of differ- 
ent states of matter, including a gas0, H 0, H- Since 
collisions with grains are inelastic, the gaseous steady 
state can only be maintained by external forcing. De- 
spite much recent work, the form of the velocity distri- 
bution for a driven granular gas remains an open ques- 
tion — velocity distribution functions found in experi- 
nwntJJlJiiJlJia IH 13 EMlja simulation 
El 111 Gj, 111 111 H3, and theory [2j H H || differ 
significantly. 

The velocity distribution function for elastic particles 
in equilibrium is Gaussian. Distributions obtained for 
inelastic granular gases are typically not Gaussian and 
are often fit to a function of the form 



P a (v) — aexp(— B\v/a\ a ) 



(1) 



where a, B and a are fitting parameters, and a = (v 2 ) 1 / 2 
H El El [H E E |U . Several experiments and simu- 
lations with different geometries and forcing mechanisms 
have found a ss 1.5 0,0,0,0] although Gaussian dis- 
tributions (a — 2) have also been observed H0. Simu- 
lations have revealed distributions that are not described 
by a single function but instead display a crossover from 
a = 2 at low velocity to a < 2 at high velocities [iH l2fl| . 
The value a = 3/2 has been obtained for the large veloc- 
ity limit for the special case of a gas of inelastic friction- 
less particles with homogeneous stochastic forcing and no 
gravity [2l| . 

Many experiments have been conducted on monolayers 

of particles [IHIlElEilElElEl 

because limiting the 
motion in one dimension allows the use of a video camera 
to record the entire velocity field. Since velocities in the 
suppressed direction can never be fully eliminated, these 
systems are quasi-2-dimensional (2D), not strictly 2D. In 
such confined geometries, particles can make as many or 
more collisions with the wall as with other grains dur- 
ing one driving cycle. Collisions with walls may then 
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influence the shape of the velocity distribution function. 
We find that the confining sidewalls play a major role 
in determining the velocity distribution function, which 
we obtain from experiments and simulations on a ver- 
tically oscillating monolayer of spheres whose motion is 
suppressed in one horizontal direction. 

Experiment. Our experimental setup, which is similar 
to that in E3, used N = 130 stainless steel balls of di- 
ameter d — 1.6 mm, contained between vertical sidewalls 
(Plexiglas plates) with a separation of l.ld. The con- 
tainer had an interior horizontal dimension of 48c? and 
vertical dimension of 32gL It oscillated with a frequency 
/ = 50 Hz, and the peak non-dimensional acceleration 
was r = Air 2 f 2 A/ g — 20, which corresponds to an am- 
plitude A = 1.25d. The container was evacuated (8 Pa) 
to avoid hydrodynamic interactions. Each run used new 
balls that were cleaned in ethanol and sonicated. In our 
experiment particles gained energy only through colli- 
sions with the bottom of the container (in earlier mea- 
surements at T = 50 EH, particles collided with the top 
as well as the bottom of the container, which was not 
evacuated). Particle motions were recorded by a digital 
camera (Phantom v4, Vision Research) at a rate of 1000 
frames per second. Particle displacements were resolved 
with an accuracy of 0.004 mm (0.0025d). Statistical 
properties were obtained by averaging over 7650 drive 
cycles and 20 different phases in the cycle. 

Simulation. An event-driven algorithm described in 
|25| was used for the simulation, which was conducted 
for the same T, /, and sidewall separation as the exper- 
iment. The parameters characterizing ball-ball interac- 
tions were the minimum coefficient of restitution e = 0.7, 
the coefficient of sliding friction /it> = 0.5, and the rota- 
tional coefficient of restitution /3 = 0.35. The coefficient 
of restitution varies with relative normal velocity (v n ) as 
described in J25| : the coefficient of restitution is the max- 
imum of e and 1 - (1 - e){vJ^gdfl A . The TC model 
of Luding and McNamara [2|j was also used to prevent 
inelastic collapse by setting the coefficient of restitution 
to unity if a particle was involved in another collision 
within 3.7 x 10 -4 seconds of the previous one. For inter- 
actions between balls and the container (both the side- 
walls and bottom), we used the same values for e and 
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/3, but we varied the coefficient of sliding friction with 
the wall from fj, w = (no sidewall or bottom friction) to 
fi w = 1. To reproduce the experiment, N = 130 parti- 
cles were simulated in a box of height 200<i, width 48d, 
and plate separation l.ld. The entire box was oscillated 
vertically so the particles collided with moving sidewalls, 
in addition to the bottom, as in the experiment. 

Steady-state distributions. Collisions of particles with 
the bottom plate inject energy mainly into vertical mo- 
tion. Energy is transferred into the horizontal direction 
directly through particle-particle collisions and through 
collisions of rotating particles with the bottom. Close to 
the bottom plate the areal density p and the probability 
distributions for the horizontal and vertical components 
of velocity (v x and v z ) vary considerably during each os- 
cillation of the plate. However, far above the plate the 
density and velocity distributions become time indepen- 




FIG. 1: (a) Linear and (b) logarithmic plots of the veloc- 
ity distribution P(v x /a x ) measured in the steady state region 
for a system with N = 130, / = 50Hz and T = 20. Distri- 
butions are shown for clean particles (•) and particles with 
a small amount of added graphite (*). Also shown is the 
non-Gaussian result (eq. 1, dashed line) from the experiment 
in [l5|. with a = 1.51 and B — 0.8. (c) Comparison be- 
tween experiment (•) with particles and simulation (o) with 
fi w = 0.075. (d) Comparison of simulations with ball-ball fric- 
tion pi, = 0.5 and with ball-wall friction y, w — 0.075 (o) and 
Pw = ("&). The experimental distributions are not precisely 
symmetric about v x = due to the container tilting slightly 
when shaking. To match the asymmetry in the experiment, 
gravity in the simulation was tilted 1.9 degrees with respect 
to the normal to the top of the container. This does not af- 
fect the functional form of the distributions when compared 
to simulations without the tilt. 



dent, as has been shown by Moon et al. |20|. Here we 
examine distribution functions for lid < z < 12d, which 
is in the steady state region - the density and horizontal 
velocity distribution functions change by less than 5% 
during each cycle. 

Our measured and simulated distributions are shown 
in Fig. ^ For clean particles in the experiment and non- 
zero wall friction in the simulation, the velocity distri- 
butions have an unusual characteristic: a sharply peaked 
maximum, a feature that has been observed before [7llll| 
but has not been fully explored. (For z < lid, the shape 
of the distribution changes slightly with height in the box 
and phase of the driving cycle, but the sharp peak is al- 
ways present.) We find that the peak disappears when 
we add approximately 0.0002 kg of gra phit e powder (a 
lubricant) to the 0.1 kg of steel spheres [22]]. The distri- 
butions observed with and without graphite both differ 
from those in [H| (cf. Fig. []]»). 

Experiment and simulation are compared in Fig. ^c). 
For fj, w — 0.075, the simulation results agree well with 
the experiment. The peak of the velocity distribution in 
the simulation decreases as p w is decreased, and the peak 
disappears completely for fi w — 0, as Fig. dd) shows. 
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FIG. 2: Double log plot of the velocity distribution functions. 
P a (v) has slope —a; to guide the eye, slope a = 3/2 is shown 
by a dashed line and a = 2 (a Gaussian) by a solid line. 
Experimental results are shown for clean stainless steel (SS) 
particles (•) and SS particles with graphite added (*). Simu- 
lation results are shown for three different values of ball-wall 
friction fi w with the ball-ball friction held fixed, /if, = 0.5. 
The data sets have been offset for clarity. 

The distributions obtained from experiments on stain- 
less steel particles with graphite and simulations with 
fi w = are described by a straight line on graphs like 
those in Fig. 2. The slope of such a graph yields the 
magnitude of the exponent a in P a (v). In simulations 
without sidewall friction, ji w = 0, the exponent obtained 
is 1.8. An exponent of 1.7 is found for the velocity dis- 
tribution of stainless steel particles with graphite. The 
peaked distributions are not described by a single value 
of a, but we can compare estimates of a local value of 
a in the range 1.0 < ln(v x /a x ) < 1.6: we obtain 1.8 
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for clean stainless steel particles, while in the simulation, 
a increases from 1.3 with sidewall friction ji w — 0.2 to 
a = 1.8 with [i w =0.1. 




FIG. 3: Horizontal velocity of a single ball on a vertically 
oscillating plate in (a) simulation (fj, w = 0.4 and fib = 0.5) 
and (b) experiment. The larger, less frequent jumps are the 
result of collisions with the bottom plate; these collisions are 
indicated by the dotted vertical lines. The more frequent 
smaller changes are the result of collisions with the sidewalls; 
these changes are clearer in the enlarged scale of the insets. 
The regions depicted in the insets are indicated by arrows. 

Single Particle Dynamics. We have shown that inter- 
actions with sidewalls strongly affect the functional form 
of the velocity distribution. This result is supported by 
our observations in simulations that in the steady state 
region, a ball collides with the wall typically three times 
as often as it collides with another ball. To isolate the 
effects of ball-wall collisions, we have conducted exper- 
iments and MD simulations on a single particle in an 
oscillating container (Fig. [^J. Because there are no col- 
lisions with other particles, the particle's motion is de- 
termined only by collisions with the bottom plate and 
the sidewalls. Figure Ola) shows the time evolution of 
the horizontal velocity v x for a particle in a simulation 
with fi w — 0.4 and e = 0.7. Each time a particle bounces 
on the bottom plate, some of the angular momentum 
of the particle can be transferred into linear momentum 
in the horizontal direction. These collisions would pro- 
duce the only changes in v x if there were no interaction 
with the sidewalls, but Fig. OJa) reveals more frequent 
smaller changes, which correspond to collisions with the 
sidewall. The staircase-like decrease in velocity (see in- 
set) corresponds to a particle's rattling between the side- 
walls, losing energy at every collision. Thus the effect of 
the sidewalls is to damp the horizontal velocity. The ul- 
timate fate of a single particle, regardless of its initial v x , 
is to bounce vertically on the bottom plate with v x = 0. 

The horizontal velocity v x measured for a single par- 
ticle in the experiment is shown in Fig. &b). Collisions 
with the bottom plate, determined to be when the verti- 
cal component of velocity v z changes sign, are indicated 



by the dotted vertical lines. If there were no influence 
of the sidewalls, the horizontal velocity v x would remain 
constant between these lines. The behavior of the par- 
ticle between collisions with the bottom plate is more 
complicated than in the simulation, but it is still clear 
that the horizontal velocity is damped by collisions with 
the walls. The damping of the horizontal motion of a 
single particle illustrated by Fig. [3] explains why the ve- 
locity distribution for a gas of particles has a peak at 
v x = (Fig. The over-populated high energy tails 
arise because for a distribution with a given variance, 
the increase in the central peak must be balanced by an 
increase for v > a. 




FIG. 4: Velocity distribution P(v/a) for a model of a damped 
driven particle. The distributions are shown for increasing 
damping 7: 0, 0.005, and 0.010. The inset compares the tails 
of the distributions for the model with a Gaussian distribution 
(solid line, a — 2) and a distribution with a = 1.5 (dashed 
line). The data sets are offset for clarity. 

Single Particle Model. Features of the velocity dis- 
tributions obtained from experiment and simulation are 
well described by a discrete map model with a damped 
driven single particle. The particle's velocity is initially 
drawn from a Gaussian distribution of variance unity. 
The velocity at iteration n+1 is given by v„+i = f n e -7 . 
For one percent of the iterations, randomly selected, we 
replace the velocity with a velocity drawn from a 

Gaussian distribution with variance unity. The veloc- 
ity probability distribution is constructed from n = 10 
iterations. The exponential decay of the particle veloc- 
ity between iteration steps corresponds to the numerous 
sidewall collisions that occur between excitations by the 
plate, and the random replacements of the particle's ve- 
locity mimic plate collisions that transfer horizontal mo- 
mentum to the particle. 

This model captures the qualitative behavior of the 
velocity distributions found in both experiment and sim- 
ulations, as Fig. 01 illustrates. For finite damping, 7 > 0, 
the distribution is strongly peaked at v = 0, while in the 
absence of damping, 7 = 0, the distribution is Gaussian. 
Further, damping affects the tails of the distribution: as 
damping is decreased to zero, double logarithmic plots 
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of the distribution become less curved and the slope in- 
creases from 1.3 to 2, just as in the MD simulation (Fig. 
0). 

The single particle model is similar to a model by 
Puglisi et al. 24] that includes damping of the parti- 
cle velocities. Increasing the damping in their model also 
led to non-Gaussian velocity distributions, but a strong 
peak around v = was not reported. This peak might 
be absent in their model because particles were driven 
not by discrete heating events but by continuous white 
noise, which for strong damping led to Gaussian behavior 
around v = in their model. 

Conclusions. The kinetic theory of granular gases is 
often studied in experiments on confined monolayers of 
grains because the behavior of all grains for all times 
can be recorded. However, we have found that the 
ball-wall friction associated with the confinement should 
be included in interpreting experiments on monolayers 
in quasi-2D geometries, including vertical jl5| . inclined 
[Tl| . and horizontal layers 0, 0. Indeed, in an exper- 



iment with the last geometry the velocity distribution 
was peaked for a smooth plate 0, but the peak dis- 
appeared when the smooth plate was replaced with a 
rough plate, which drove horizontal as well as vertical 
motion |9j . Similarly, a recent experiment with a layer of 
light particles on top of a layer of heavy particles yielded 
a non-Gaussian distribution for the heavier particles, but 
Gaussian statistics were found for the lighter particles |(| . 
The interactions between the particles and the container 
in these quasi-2D systems may have been principal deter- 
minants of the shape of velocity distributions and there- 
fore should be taken into consideration. 
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